function [mutually_acceptable] = get_consideration_sets(couple,singlemale,singlefemale,ageM,ageF,lb,ub)

n = length(couple);
mutually_acceptable = zeros(n,n);

for i=1:n
    if couple(i) == 1
        mutually_acceptable(i,i)=1;
    end
    for j=1:n
        if ((couple(i) == 1 && couple(j) == 1) && (ageM(i)- ageF(j) <= ub) && (ageM(i) - ageF(j) >= lb))
            mutually_acceptable(i,j)=1;
        end
        if ((couple(i) == 1 && singlefemale(j) == 1) && (ageM(i)- ageF(j) <= ub) && (ageM(i) - ageF(j) >= lb))
            mutually_acceptable(i,j)=1;
        end
        if ((singlemale(i) == 1 && couple(j) == 1) && (ageM(i)- ageF(j) <= ub) && (ageM(i) - ageF(j) >= lb))
            mutually_acceptable(i,j)=1;
        end        
     end
end

return